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Abstract 

Two species of animals are competing in the same environment. 
Under what conditions do they coexist peacefully? Or under what 
conditions does either one of the two species become extinct, that is, 
is either one of the two species excluded by the other? It is natural to 
say that they can coexist peacefully if their rates of reproduction and 
self-limitation are relatively larger than those of competition rates. In 
other words, they can survive if they interact strongly among them- 
selves and weakly with others. We investigate this phenomena in 
mathematical point of view. 

In this paper, we concentrate on coexistence solutions of the compe- 
tition model 

Au + u{a — g{u,v)) = . ^ 
< Av + v{d- h{u,v)) =0 ™ ' 
u\dn = v\d^i = 0, 

This system is the general model for the steady state of a competi- 
tive interacting system. The techniques used in this paper are ellip- 
tic theory, super-sub solutions, maximum principles, implicit function 
theorem and spectrum estimates. The arguments also rely on some 
detailed properties for the solution of logistic equations. 



1 Introduction 

A lot of research has been focused on reaction-diffusion equations model- 
ing of various systems in mathematical biology, especially the elliptic steady 
states of competitive and predator-prey interacting processes with various 
boundary conditions. In the earlier literature, investigations into mathemat- 
ical biology models were concerned with studying those with homogeneous 
Neumann boundary conditions. From here on, the more important Dirichlet 
problems, which allow flux across the boundary, became the subject of study. 

Suppose two species of animals, rabbits and squirrels for instance, are com- 
peting in a bounded domain Q. Let u{x, t) and v{x^ t) be densities of the two 
habitats in the place a; of ^2 at time t. Then we have the following biological 
interpretation of terms. 

[A) The partial derivatives Ut{x,t) and Vt{x,t) mean the rate of change of 
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densities with respect to time t. 

(B) The laplacians Au{x,t) and Av{x,t) imply the diffusion or migration 
rates. 

(C) The rates of self-reproduction of each species of animals are expressed as 
multiples of some positive constants a, d and current densities u{x^ t),v{x, t), 
i.e. au{x, t) and dv{x, t) which will increase the rate of change of densities in 
{A), where a > 0,d > are called self-reproduction constants. 

(D) The rates of self-Hmitation of each species of animals are multiples of 
some positive constants b, f and the frequency of encounters among them- 
selves u'^{x,t),v'^{x,t), i.e. bu'^{x,t) and v'^{x,t) which will decrease the rate 
of change of densities in (^4), where 6 > 0, / > are called self-limitation 
constants. 

(E) The rates of competition of each species of animals are multiples of some 
positive constants c, e and the frequency of encounters of each species with 
the other u{x,t)v{x,t), i.e. cu{x,t)v{x,t) and eu{x,t)v{x,t) which will de- 
crease the rate of change of densities in (^4), where c > 0, e > are called 
competition constants. 

(F) We assume that both species of animals are not staying on the boundary 
of Q. 



Combining all those together, we have the following dynamic model 

Ut{x, t) = Au{x, t) + au{x, t) — bu'^{x, t) — cu{x, t)v{x, t) 
Vt{x, t) = Av{x, t) + dv{x, t) — ] 
u{x, t) = v{x, t) = for x e dD,, 

or equivalently, 



Vt{x, t) = Av{x, t) + dv{x, t) — fv'^{x, t) — eu{x, t)v{x, t) ^ ^ °°)) 



Ut{x,t) — Au{x,t) + u{x,t){a — bu{x,t) — cv{x,t)) • o rn ^ 
Vt{x,t)=Av(x,t)+v{x,t)(d-fv{x,t)-eu{x,t)) i^x[U,ooj, 
u{x,t) — v{x,t) = for x e dfl, 

Here, we are interested in the time independent, positive solutions, i.e. the 
positive solutions u{x),v{x) of 

Au{x) + u{x){a — bu{x) — cv{x)) = • q 
Av{x) + v{x){d - fv{x) - eu{x)) =0 ' 
. Mdn = v\oa = 0, 
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which are called coexistence state or steady state. The coexistence state is 
the positive density solution depending only on spatial variable not on 
time variable t, and so the existence of that means the two species of animals 
can live peacefully and forever. 

A lot of work about the existence and uniqueness of coexistence state of 
the above steady state model has already been established during the last 
decade. (see [1], [2], [3], [6], [7], [9], [10].) 

In this paper, we study rather general types of system. We concern the 
existence and uniqueness of positive coexistence when the relative growth 
rates arc nonlinear, more precisely, the existence and uniqueness of positive 

steady state of 

Au + u{a — g{u,v)) — . „ 
< Av + v{d-h{u,v))^0 ' 
u\dn = v\dn = 0, 

where a, d are positive constants, g, h are functions, is a bounded do- 
main in i?" and u, v are densities of two competitive species. 

The foUowings are questions raised in the general model with nonlinear 
growth rates. 

Problem 1 : Under what conditions do they coexist? Under what conditions 
do they have unique steady state? When does either one of the species be- 
come extinct? 

Problem 2 : Assuming that they can coexist and the coexistence state is 
unique at a fixed self- reproduction (a, d), can they still coexist regardless of 
slight change of that self-reproduction? 

Problem 3 : This is a generalization of Problem 2. If we have the existence 
and uniqueness of coexistence state on the left boundary of a closed convex 
region F for the reproduction (a, d), can we extend the region F to an open 
set including F without losing the uniqueness? 



In section 3, some sufficient conditions to guarantee the existence and unique- 
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ness of positive solutions were obtained, and we could also see that there is no 
positive solution for small self-reproduction rates mainly using upper-lower 
solutions and spectrum estimates, which solves Problem 1. In sections 4 and 
5, we answer problems 2 and 3 using elliptic theory, maximum principles and 
implicit function theorem. 

2 Preliminaries 

In this section we will state some mathematical preliminary results which 
will be useful for our later arguments. 

Definition 2.1 (Super and sub solutions) The vector functions {u^, ...,u^) 
and {y}, form an super/sub solution pair for the system 

Au' + g\u^, ...,u^) = m n 
= on dfl 

if fori ^ 1,...,N 

' Au' + g'{u\ M^-i, u\ u'+^, ...,u^) <0 
< Au' + g'{u\...,u^-\u\u'^+\...,u^) > 
in fl for yP < < u^,j ^ i, 

and 

< -u* on Q 

< < -u* on d^l. 

Lemma 2.1 If g^ in the Definition 2.1 are in and the system admits an 
super/sub solution pair {u^, ...,u^), {u^, ...,u^), then there is a solution of 
the system in 2.1 with < ?/ < in Cl. If 

Au' + g\u\...,u^) ^ 0, 
Au' + g\u^,...,u^) 7^ 

in Q for i = 1, N , then < < in 

Lemma 2.2 (The first eigenvalue) 

—Au + q{x)u = Xu in fl, , . 

u\9a = 0, ^ ^ 
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where q{x) is a smooth function from fl to R and fl is a bounded domain in 

(^4) The first eigenvalue Xi{q) of (1), denoted by simply Ai when q = 0, is 

simple with a positive eigenfunction. 

{B) If qi{x) < q2{x) for all x G il, then Ai(gi) < Ai(g2)- 

{C) (Variational Characterization of the first eigenvalue) 

\i[q) = mm 



Lemma 2.3 



= ^ aij{x)DijU + '^ai{x)DiU + a{x)u = f{x) in ft, 

i,j=l 1=1 

where Q is a bounded domain in i?" and 
(Ml) dn G C2'°(0 < a < 1), 

(M2) \aij{x)\a, \ai{x)\a, \a{x)\a <_M{i,j = 1, ...,n), 

(M3) L is uniformly elliptic in VL, with ellipticity constant 7, i.e., for every 
X G and every real vector ^ — (^1, ^n) 

n n 
i,j=l i=l 

1. Maximum principles 

Let u G C^{Q) n C{Q) be a solution of Lu > 0{Lu < 0) in fl. 

(A) If a{x) = 0, then maxf^ u — max^o u{mmQ u — min^n u) . 

(B) If a{x) < 0, then maxQU < maxgQU'^{mmQU > —maxgciW), 

where = max(-u, 0),u~ = — min(-u, 0) . 

(C) If a{x) = and u attains its maximum (minimum) at an interior point 
offl, then u is identically a constant in Q. 

(D) If a{x) < and u attains a nonnegative maximum (nonpositive mini- 
mum) at an interior point offl, then u is identically a constant in fl. 

2. Schauder's estimate 

Ifue C^'^'ifl) andu\aa^(l)e C^'°'{dfl), then 



U\2,a < c{\Lu\a + \u\o + \<r\2,aJ: 



where the constant c> is independent of u. 
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Lemma 2.4 (Implicit Function Theorem) 

Let X, Y, Z be Banach spaces. For a given (mq, vq) & X xY and a,b > 0, let 
S — {{u,v) :|| u — uo \\< a,\\ V — vq \\< b}. Suppose F : S ^ Z satisfies the 

following: 

{A) F is continuous. 

[B) Fy{-, ■) exists and is continuous in S(in the operator norm). 



{D) [Fy{uQ,vo)]~^ exists and is a continuous map from Z to Y. 
Then there are neighborhoods U of uq and V of vq such that the equation 
F{u,v) = has exactly one solution v &V for every u & U and the solution 
V depends continuously on u. 

We also need some information on the solutions of the following logistic 
equations. 

Lemma 2.5 (in [9]) 



where Q is a bounded domain in and 

{A) f is a strictly decreasing function, 

(B) there exists Cq > such that f{u) < for u > Cq. 

(1) ///(O) > Ai, where Ai is the first eigenvalue of —A with homogeneous 
boundary condition, then the above equation has a unique positive solution. 

(2) ///(O) < Ai, then u = is the only nonnegative solution of the above 
equation. 

In the case (1), we denote this unique positive solution as 6f. The main 
property about this positive solution is that 6'/ is larger as / is larger, i.e. 
0g<9fiig< /.(Proposition ??) 



(C) F{uo,Vo) = 0. 




An + uf(u) =0 m n, 
u\dii = 0,M > 0, 
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3 Existence, Nonexistence and Uniqueness of 
steady state 

We consider the elliptic system 

Au + u{a — g{u,v)) = . „ 

Av + v{d-h{u,v))^0 ' (2) 
u\dn = v\9^i = 0. 

Here Q is a bounded, smooth domain in i?" and 

(Ul) g,h E arc strictly increasing functions with respect to u,v, 
{U2) there exist ki, k2 > such that g{u, 0) > a for u > ki and 
h{0, v) > d for V > /c2- 

If there is no competition between the species, that is, if we consider 

Au + u{a- g{u,0)) = . 
Av + v{d-h{0,v)) = ' 
u = V = on dQ, 

then by the Lemma 2.5, the condition a > \i,d > Ai(i.e. reproductions are 
relatively large.) were sufficient to guarantee the positive density solution 
9a-g{-fi),9d-h{o,-)- But, if there is some competition between them, then as 
see in the following Theorem 3.1, we should have larger lower bound for repro- 
duction rates a and d, i.e. we have stronger conditions a > Ai + g{0, k2),d > 
Ai + h{ki, 0) to guarantee their coexistence. (i.e. the reproductions should be 
much larger.) 

The following is the main result: 

Theorem 3.1 (A) If a > Ai + g{0, k2) and d > Ai + h{ki, 0), then (2) has a 
positive solution {u, v) with 

Conversely, any positive solution {u,v) to (2) must satisfy the inequalities. 
(B) //a > Ai + 51(0, k2) and d> \i + h{ki, 0) and 

^^^n^ inf — > sup sup — + sup sup — 

.dg. ,dh. 
+ 2sup(-)sup(-), 
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where B = [0, ki] x [0, k2], then (2) has a unique coexistence state. 
iC) If a < Xi or d < Xi, then (2) does not have any positive solution. 

Biologically, we can interpret the conditions in Theorem 3.1 as follows. The 
constants a, d and functions (?, h describe how species 1 {u) and 2 {v) interact 
among themselves and with each other. Hence, the both conditions in [A) 
and {B) imply that species 1 interacts strongly among themselves and weakly 
with species 2. Similarly for species 2, they interact more strongly among 
themselves than they do with species 1. The inequalities in the conclusion 
{A) imply that the densities with competitions(-u and v) are less than those 
without competition. (6'a_(,(.^o) and 6*^-/1(0,-)) Furthermore, (C) says that if one 
of the species has small reproduction, then it may be extinct, which means 
that the two species can not coexist. 

Proof. [A) Let u — 9a-g(.,o),v — Od-h{o,-)- Then since g is increasing, we 



have 



Au + u{a — g{u, v)) 

Am + u{a — g{u, 0) + g{u, 0) 

uigi'^: 0) - gi^: ^)) < 0- 



9{u,v)) 



Similarly, we have 

Av + v{d- h(u,v)) < 



So, {u,v) is a super solution of (2). 

Let U = 0a-g{-,k2)^^^ Od-h{ki,-)- 

Then by the Maximum Principle, we obtain 




Since g is increasing, we get 



Au + u{a 
Au + u{a 
u{g{u, /C2) 



g{u, /cs) + g{u, /ca) 
g{u,v))>0. 



g{u, V)) 



Similarly, we get 

Av + v{d - h{u,v)) > 0. 



Therefore, {u,v) is a lower solution of (2). 

Furthermore, u<u and v<vmQ and u = u = v = v = Oon dfl. 
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So, by Lemma 2.1, (2) has a solution {u,v) with 

Suppose {u,v) is a positive solution to (2). By the Mean Value Theorem, 
there is v* such that 

g[u, v) = g[u, 0) H — v. 

Then 

doiu v*^ 

Au + u{a - g{u, 0)) = — -uv > in Q. 

Hence, m is a subsolution to 

f Az + z{a-g{z,0))^0 in Q, 
[ z\dn = 0. 

Any sufficiently large positive constant is a super solution to 

f Az + z{a- g{z,0)) = in Q, 
\ z\dn = 0. 

Therefore, by the super-sub solution method, we have 

U < Oa-g(.,0)- (3) 

The same argument shows 

V < Od-h{o,-)- (4) 

For sufficiently small e > 0, 

e^6'd-/i(o,-) + <^Od-h{o,-)id — h{0, e9d-h{o,-))) 
= €[A9d-h(o,-) + ^d-h(o,-){d - /i(0, e9d-h(o,-)))] 
> e[A9d-h{o,.) + 9d-h{o,-){d - h{0, 9d-h{o,-)))] 
= in Q, 

and so e9d-h{o,-) is a sub solution to 

f Az + z{d- h{0,z)) = in fl, 
\ z\q^ = 0. 

10 



Since d — h{0, k2) < 0,k2 is a super solution to 

f Az + z{d-h{0,z))^0 in 
\ z\on = 0. 

Hence, by the super-sub solution method again, 
So, 

g{u,v) < g{u,0d-h{o,-)) < g{u,k2) 

since g{u,z) is increasing. Therefore, 

Au + u{a — g{u, k2)) < Au + u{a — g{u,v)) 

= in 

Hence, w is a super solution to 

{Az + z{a — g{z,k2)) = in Q, 
z\dn = 0. 

Let 4>i be the first eigenvector of 

f Am + Aiz = in Q, 
\ ^|af2 = 0. 

Then for sufficiently small e > 0, 

a — g{e4>i, k2) — Ai > in O, 

and 

A(e0i) + e0i(a - 5((e0i, A;2)) = e[A(e0i) + 0i(a - ^((e^i, ^2))] 

e(A0i + Ai0i) = in Q. 

Consequently, e0i is a sub solution to 

J Az + 2;(a — 51(2;, A;2)) = in Q, 
I z\dQ, = 0. 
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Hence, by the super-sub solution method again, 

da-gi-M) < (5) 

The same argument shows 

Od-h(ki,-} < V. (6) 

Prom (3) to (6), we have 

da-gi-M) < Oa-g(.fi), Od~h{ki,-) < V < Od-h{0,-) (7) 

Consequently, for any positive solution {u,v) of (2), the inequalities (7) hold. 

(B) Suppose {ui,vi) and {u2,V2) are positive solutions to (2). 
Let p = Ui — U2 and q = Vi — V2. Then 

Ap+ {a- g{ui,vi))p = Aui - Au2 + {a - g{ui,vi)){ui - U2) 

= -Au2- {a- g{ui,vi))u2 
= -Au2 - U2{a - g{u2, V2) + g{u2, V2) - g{ui, vi)) 

= -U2{9{U2:V2) - giUi,Vi)) 

= -U2{g(u2, V2) - g{ui, V2) + g{ui, V2) - g{ui, vi)) 
dg{x,V2) , . dg{ui,x) 

= + inn, 

where from Mean Value Theorem depending on ui, U2, Vi,V2. Hence, 

Ap + (a- v,))p - u,lj,?^ + 5^^) = in (8) 
Similarly, we can get 

Aq + {d-h{u2,V2))q-Vi{p — \- q — ) = m Q, (9) 

where y,y are from Mean Value Theorem depending on Ui,U2,Vi,V2. Since 
Ai(a — g{ui,Vi)) = 0, by the Variational Characterization of the first eigen- 
value, 

z{-Az - {a- g{ui,vi))z)dx >0 (10) 



/ 

Jn 
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for any z G C^(f2) and z\q^ = 0. The same argument shows that 

j^w{-Aw - {d- h{u2,V2))w)dx>0 (11) 

for any w e C^(fi) and w\qq — 0. Prom (8) and (9), we have 

r _p^p _ _ g(^u„ v,))p^ + U2p{p'^^ + = 

1 -gAg - (rf - ^2))g^ + v,q{p'-^ + q'-^) = 

Using (10) and (11), we have 

f dg{x,V2) dg{ui,x) dh{y,vi) dh{u2,y).. 

Hence, 

dg{x,V2) 2 , / dg{ui,x) dh{y,vi) dh(u2,y) 21 ^ ^ 

+ (^^^^ + ^^^^)^^ + ^^^^^ 1 ^ 0- 



Therefore, p = g = if we can show that 

9^(^ ^ ^^a^^l^ _ ^^^^^ag(x,.2)9M^2,y) ^ Q 
dv ou ou dv 



which is true if 



I.e., 



or 



-Wi^^^^<Oina 



19., „ dh{y,vi) ■ Q 

^ ag(a,t)2) dh{u2,y) ^ M2 / 8g("i,^) \2 , fi / '\2 
du dv vi^ dv ^ U2^ du ' 



I 2 d9{ui,x) dh(y,vi) q 
dv du ' 

This is the case from the hypothesis in the theorem and (7), and so the 
uniqueness is proved. 
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(C) Assume a < Ai. Suppose {u^v) is a nonncgative solution to (2). Then 
since g is an increasing function with respect to u and v, 

Am + u{a — g{u, 0)) 
= Au + u{a — g{u,v) + g{u,v) — g{u,0)) 
= u{g{u,v) - g{u,0)) >0. 

Therefore, m is a sub solution to 

J Au + u{a — g{u,Q)) = Q in 
\ u\dn = 0. 

Any constant larger than ki is a super sohition to 

f Au + u{a-g{u,0)) ^0 in 
\ -"Ian = 0. 

Hence, by the Lemma 2.1, there is a solution u of 

f Aii + ii(a - g(u, 0)) = in Q, 
\ u\dn = 

such that < u < u. But, since a < Ai, = by (2) of Lemma 2.5, and so 
u = 0. 



4 Uniqueness with small perturbation of re- 
production rates 

We consider the model 

Au + u{a — g{u,v)) = . 

Av + v{d-h{u,v)) = ' (12) 
u\dn = vlon = 0. 

Here f2 is a bounded, smooth domain in i?" and 

(PI) g,h E are strictly increasing functions with respect to u and v, and 
^(0,0) = M0,0)=0, 

(P2) there are ki,k2 > such that g{u, 0) > a > Xi ior u > ki and /i(0, v) > 
d > Xi ior V > k2- 

The following is the main theorem. 
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Theorem 4.1 Suppose 

{A) a > \i{g{0, Od-hio,-))), d > Xi{h{9a-g(.,o), 0)), 

(B) (12) has a unique coexistence state {u,v), 

(C) the Frechet derivative of (12) at {u,v) is invertible. 

Then there is a neighborhood V of (a, d) in such that if (oq, do) G V, then 
(12) with {a,d) = (ao,(io) has a unique coexistence state. 

Theorem 4.1 looks hke the consequence of Imphcit Function Theorem. But 
the inverse function theorem only guaranteed the uniqueness locally. Theo- 
rem 4.1 concluded the global uniqueness. The techniques we will use includes 
naturally Implicit Function Theorem and a priori estimates on solutions of 
(12). 

Biologically, the first condition in this theorem indicates that the rates of 
self-reproduction is large. The condition of invertibihty of Frechet derivative 
also illustrates that the rates of self-limitation is relatively larger than those 
of competitions which will be in Theorem 4.3. Then the conclusion says 
that small perturbation of reproduction rates does not lose the existence and 
uniqueness of positive steady state, i.e. they can still coexist peacefully even 
if there is some slight change of reproduction rates. 

Proof. Since the Frechet derivative of (12) at {u,v) is invertible, by the 
Implicit Function Theorem, there is a neighborhood V oi (a, d) in and a 
neighborhood W of {u,v) in [Cq^"" {fl)]'^ such that for aU (oq, (io) e V, there 
is a unique positive solution {uq,vq) E W of (12). Suppose the conclusion of 
the theorem is false. Then there are sequences (a„, dn, Un, Vn), {an, dn, u^, v*) 
in y X [Cq'^°'{Ci)]'^ such that (m„,v„) and (M*,f*) are the positive solutions 
with {a,d) = {an,dn) and {un,Vn) ^ «,<) and (a„,d„) {a,d). By the 
standard elliptic theory, (ti„,f„) {u,v) and {u^,v!^) {u*,v*) in C^'", 
and {u,v), {u*,v*) are solutions of (12). Claim -u > 0, -u > 0, -u* > 0, f * > 0. 
It is enough to show that u and v are not identically zero because of the 
Maximum Principle. Suppose not, then by the Maximum Principle again, 
one of the following cases should hold: (1) uis identically zero and v > 0. (2) 
u > and v is identically zero. (3) u is identically zero and v is identically 
zero. 

Without loss of generality, assume u is identically zero. 
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Let Un = II "Yi , Vn = Vn for all n E N. Then 



\ AVn + Vn{dn " h{Un, Vn)) = 

From the elliptic theory, Un —>■ u and 

r Ati + t}^(a-(?(0,t;)) = 

I Ai; + ^;(d-/i(0,^;)) = ' 

since g,h are continuous, i.e., a = \i{g{0,v)). 

(1) If V = 0, then by the monotonicity of g and Ai, a = Xi[g{0,v)) — 
\i{g{0,0)) < Xi{g{0,9d-h{o,-))) which contradicts our assumption. 

(2) If V is not identically zero, then v = 9d-h{o,-) and so a = Xi{g{0,v)) = 
Xi{g{0,9d-h{o,-))) which is also a contradiction to our assumption. Conse- 
quently, {u,v) and {u*,v*) are coexistence states for {a,d). But, since the 
coexistence state with respect to {a,d) is unique, {u,v) = {u*,v*) = {u,v). 
But, since {un,Vn) 7^ (m*,^*), it contradicts the Implicit Function Theorem. 

The proof of the theorem also tells us that if one of the species becomes 

extinct, in other word, if one is excluded by others, then that means the 
reproduction rates are small, i.e. the region condition of reproduction rates 
(A) is reasonable. 

Theorem 4.2 // (a„, dn, m„, f „) (a, d, u, v) and if u = or v = 0, then 
a < Xi{g{Q, 9d-h{o,-))) or d < Xi{h{9a-g(.,o): 0)). 

The condition, invertibility of Frechet derivative, in Theorem 4.1 is too arti- 
ficial. Now we turn out attention to get conditions to guarantee the invert- 
ibility of the Frechet derivative. 

Theorem 4.3 Suppose {u,v) is a positive solution to (12). 

//4inf ^S^inf ^^Mt; > [sup ^^^w + sup then the Frechet 

derivative of (12) at {u,v) is invertible. 

Proof. The Frechet derivative at {u, v) is 

( -A + g{u,v) + u'-^-a «^ \ 
I v'-^ -A + h{u, v) + v^ -dj- 
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We need to show that N{A) = {0} by Fredholm alternative. If 

r -Acp + (giu, v) + - a)^ + 'S^ui^ = 0, 

1 -AV' + '-^vv + {h{u, v) + v'-^ - d)i^ = 0, 



then 



SnW? + ^^V'V' + {h{u, v) + - d)r\ = 0. 

Since Ai(g'(w, v) — a) — Xi{h{u, v) — d) — 0, 



/o[|V^P + (^7Kt;)-aV2] >0, 
la[m? + {h{u,v)-d)i,']>0. 



Hence, 



InC-^W^ + '-%^vr) < 0. 

Hence, + + ^^;)v^^ + '-^v^] < 0. 

Hence, if 4 inf (^) inf(^)tzt; > [(sup(^))M + (sup(^))t;]2, then 
the integrand in the left side is positive definite form in Q, which means 
if = t/j = 0. Therefore, the above Prechet derivative A is invertible. 

Combining the Theorems 3.1, 4.1 and 4.3, we have the following which is 
actually the main result in this section. 

Corollary 4.4 Suppose 

(A) a> Xi+ g{0, k2), d> Xi + h{ki, 0), and 

{B) 

4inf^^inf^^ > [sup ^ + sup ^ sup ^f^] 
where B = [0, h] x [0, k2\. 

Then there is a neighborhood V of (a, d) in such that if {ao, do) £ V, then 
(12) with {a,d) = {ao,do) has a unique coexistence state. 
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Proof. From 6'a_c,(,o) < ki, 9d-h{o,) < A;2, and the monotonicity of 51(0, •), /t(-, 
we have 

( a> \i+ g{0, /C2) > Xi{g{0, ^d-Mo,-)))' 
\d>Xi + h{ki, 0) > Ai(/i(e„_5(,o), 0)). 

B dx B dy 



0) 



> [sup — ^ h sup — ^ — sup [sup — ^ — sup 

+ sup 



dy dx Oa-g{o,k2) dy 9d-h{ku-) 

dh{x,y) 



dx 

^ [sup ^^^f ' sup + sup ^^i^- sup 

oy Od-h(ku-) dy dx 

+ sup sup sup sup 

"d-h(ki,-) ^a~g{-,k2) dy ox 

dh{x,y)2 Od-h{o,-) 

+ [SUP ] sup- ^ 

^a-g(-,k2) 

^ r dg{x,y) 2 Oa-g(.,o) dg{x,y) dh{x,y) 

> sup — sup- — ^^ + 2 sup — sup — 

dy Od-hiki,-) dy dx 

dh{x,y)2 Od-h{o,-) 
+[sup — ] sup- ^ 

ox t^a-gi-M) 
since ^a-g(-,0) > Oa-g{-,k2)i^d-h{Q,-) > Od-h(ki,-)- 

Therefore, (12) has a unique coexistence state {u, v) from Theorem 3.1. Fur- 
thermore, by the estimate of the solution in the proof of Theorem 3.1, 

^^^^ dgix,y) .^^ dh{x,y) 
B dx B dy 

^ r dg{x,y) dh{x,y) Od-hio,-) dg{x,y) Oa-g{;0) 

> [sup — - — - + sup — ^ — - sup [sup — - — - sup ^ ' 



+ sup 



dy dx ^a-gi-M) dy Od-h{ki,-] 

dh{x,y) 
dx 



^ r dg{x,y) dh{x,y)v dg{x,y)u dh{x,y) 

> sup — h sup — sup — + sup — 

dy dx u dy V dx 
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Thus, we obtain 

cdg(x,y). dh(x,y) dg{x,y) dh{x,y) ^ 

4 mi — 7; — - mr — 7; — -uv > sup — 7; — -u + sup — 7; — -v\ . 
B ax B dy ay ox 

It implies that the Prechet derivative of (12) at {u, v) is invertible from The- 
orem 4.3. Therefore, the theorem follows from Theorem 4.1. 

5 Uniqueness in a region of reproduction rates 

Consider the model 

Am + M(a — w)) = . „ 

Av + v{d-h{u,v))^Q ' (13) 
u\au = v\dn = 0. 

Here Q is a bounded smooth domain in i?" and g,h E are strictly increas- 
ing functions with respect to u and v, and g{0, 0) = h{0, 0) = 0. 

The following is the main theorem. 
Theorem 5.1 Suppose 

(A) Fis a closed, convex region in B? such that for all (a, d) e F, 

a > Xi{g{Q,9d-hio, ))) and d > \i{h{9a-g(.fi),0)), 

(B) there exist cq > and Ci > such that for all (a, d) & T, g{x,0) > a > Ai 
for x > Cq and h{0, y) > d > \i for y > Ci, 

(C) (13) has a unique positive solution for every {a,d) G OlT, 
where diF = {(A^, d) e r\For any fixed d, A^ = inf{a|(a, d) e F}}, 

(D) for all (a, d) e F, the Frechet derivative of (13) at every positive solution 
to (13) is invertible. 

Then for all {a,d) G F, (13) has a unique positive solution. Furthermore, 
there is an open set W in such that V <^ W and for every (a, d) G W , 
(13) has a unique positive solution. 

Theorem 5.1 goes even further than Theorem 4.1 which states the uniqueness 
in the whole region of (a, d) whenever we have the uniqueness on the left 
boundary and invertibility of linearized operator at any particular solution 
inside the domain. 
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Proof. For each fixed d, let = sup{a : (a, d) G F} and = inf{a|(a, d) G 
F}. We need to show that for every a such that A^ < a < A'*, (13) has a 
unique positive solution. Since (13) with (a, d) — (A^, d) has a unique posi- 
tive solution {u,v) and the Prechet derivative of (13) at {u,v) is invertible, 
by theorem 4.1, there is an open neighborhood V of (A^^, d) in such that if 
{aQ,do) G V, then (13) with (a.d) = (ao,(io) has a unique positive solution. 
Let Xs = sup{A > Arf : (13) has a unique coexistence state for A^ < a < A}. 
We need to show that As > A'^. Suppose A^ < A*^. Prom the definition of 
Xs, there is a sequence {A„} such that A„ — * A~ and there is a sequence 
{un, Vn) of the unique positive solution of (13) with (a, d) = (A„, d). Then by 
the EUiptic theory, there is {uo,Vo) such that (M„,fn) converges to {uo,Vo) 
uniformly and (^0,^0) is the solution to (13) with {a,d) = {Xs,d). We claim 
that Uq is not identically zero and Vq is not identically zero. Suppose this 
is false. Then by the Maximum Principle, one of the following cases should 
hold: (1) Uq is identically zero and Vq is not identically zero. (2) Uq is not 
identically zero and vq is identically zero. (3) Both uq and vq arc identically 
zero. The argument is similar to what we had in the previous section. 
(1) Suppose uo is identically zero. 
Let Un — II "7, and Vn — Vn for all n E N. Then 

f AUn + M„(A„ - g{Un, Vn)) = Q 
[ AVn + Vn{d - h{Un, Vn)) = 

We know tin ^ u from the elliptic theory, and 



Au + u{Xs-g{0,vo))^0 

Avo + vo{d-h{0,vo))^0 '''''' 



since g,h are continuous. Hence, vq — 9d-h{Q,-) and A^ = Xi{g{0,vo)). If 

Vq is identically zero, then by the monotonicity of g and Ai, we have A^ = 

Xi{g{0,VQ)) = Ai (6/(0,0)) < Xi{g{0,9d-h{o,-))) < Ad which is impossible. If 

Vq is not identically zero, then vq = 6d-h{<d,-) and so A^ = X\{g{Q,VQ)) — 

Xi{g{0,9d-hio,-))) < Xd which is also impossible. 

(2) Suppose Vq is identically zero. 

Let Un — Un and Vn = . ^7 for all n G A^. Then 

W^n II oo 

f AUn + Un{Xn - g{Un, Vn)) = ^ Q 
I AVn + Vn{d - h{Un, Vn)) = 
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Again Vn ^ v hy the elliptic theory, and 

Auo + uo{Xs- g{uo,0)) = . 
/^v + v{d - h{uQ,Q)) ^ Q 

since g,h are continuous. Hence, d — Ai(/i(ito, 0)). If Uq is identically 

zero, then by the monotonicity of h and Ai, we have d = Xi{h{uQ,0)) — 
Ai(/i(0,0)) < \i{h{9x,-g(.fi),0)) < Ai(/i(^Ad-g(-,o)! 0)) which is impossible, 
since {X'^,d) G F. If Uq is not identically zero,, then Uq = 6'A^-g(.,o) and 
so, d = Ai(/i(mo,0)) = Xi{h{9x,-g(.fl),0)) < Ai(/i(6'Ad_g(.,o), 0)) which is also 
impossible, since (A'', d) e F. Consequently, uq > 0,vo > in Q, that is, 
(uo, fo) is a coexistence of (13) with (a, rf) = (As, d). Since (A^, c?) G F, by the 
assumption, the Frcchct derivative of (13) with {a,d) = {Xg.d) at (uq.Vq) is 
invertible. Hence, by the Implicit Function Theorem, there is an open neigh- 
borhood U of As and an open neighborhood V of (mq, Vq) such that if a G C/, 
then (13) has a unique coexistence state in V. But, by the definition of As, 
there is a sequence {A^} C U such that A^ — >• A+ and there is a sequence 
{{u'^,v'^)} of the coexistence state of (13) with {a,d) = (A^, ci) such that 
{u'^, v'^) ^ V for all n E N. By the Elliptic Theory again, u'^ — >• u'q, v'^ —>■ v'q 
and from the same argument above, (^0,^0) ^ V is also a coexistence of 
(13) with (a, d) = (As,d). Since (As,rf) G F, by the assumption again, the 
Frcchct derivative of (13) at {uq,Vq) is invertible. Hence, by the Implicit 
Function Theorem again, there is an open neighborhood U' of As and an 
open neighborhood V of {uq,Vq) such that if a G U', then (13) has a unique 
coexistence state in V. Consequently, there are points in the left side of As 
such that (13) has two different coexistence states. That is a contradiction 
to the definition of As. Hence, As > A'' and the first part of the theorem is 
proved. Furthermore, by the assumption, for each {a,d) G F. the Frechet 
derivative of (13) at the unique solution {u,v) is invertible. Hence, Theorem 
4.1 concludes that there is an open neighborhood V(^a,d) of {a,d) in such 
that if (oo, do) G V(^a,d), then (13) with reproduction rates {qq, do) has a unique 
coexistence state. Let W = [j{a,d)er ^ia,d)- Then W is an open set in R"^ such 
that F C and for each (oq, do) G W, (13) has a unique coexistence state. 

Apparently, Theorem 5.1 generalizes Theorem 4.1 and consequently, we have 
the following which is actually the main conclusion in this section. 
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Corollary 5.2 Suppose 

{A) V is a closed, convex region in R^, 

(B) there exist ki,k2 > such that for all {a,d) e F, a > Ai + g{0,k2), 

d> Xi + h{ki, 0), a - g{ki, 0) <0,d- h{0, ki) < 0, 

(C) 



A\-nf Mflllinf 9h{x,y) r dg{x,y) , dh{x,y) 6 

4 mts mtB —Q— > [sup —g— + sup sup(„ ,^)gr j 



d-h(0,-) 



dy " i^^f dy ' dx •^"i^(M)er e,_,(.,fc2)J 

[sup^sup(„,,),,^ + sup^], 

where B = [0, /ci] x [0, k2\ ■ 

Then there is an open set W in E? such that T C.W and for every (a, d) e V, 
(13) has a unique positive solution. 

The condition {B) means V is some set of large self-reproduction rates, and 
the condition (C) implies that the self-limitation rates are relatively larger 
than competition rates. Then the conclusion says that the existence and 
uniqueness of coexistence state are guaranteed on F and the region F can be 
extended to a larger set without losing the uniqueness. 

Proof. From 9a-g{.fl) < ki, Od-h{o,-) < ^2, and the monotonicity of 
g{0, •), 0) we have 

r a > Ai + ^(0, /C2) > Ai(^(0, 9d-h{o,.))), 
\ d > Ai + h{ki, 0) > Ai(/i(^„_^(.,o), 0)). 

for all (a, d) e F. 

By the condition (C), for every (a, d) e 9F, 

4inf M^inf 

B dx B dy 

^ r dg{x,y) dh{x,y) ^d-Mo,-) u dg{x,y) ^a-a(-,o) 
> [sup — - — - + sup — - — - sup ^-^J [sup — - — - sup 



-|- sup 



dy dx Oa-g(o,k2) dy Od-h{ku-) 

dh{x,y) 
dx 



r dg{x,y) 2 ^a-s(-,o) . dg{x,y) dh{x,y) 
sup — sup- — ^^+sup — sup — 

dy Od-h(ki,-) dy dx 

0a-g{-fi) dd~h(o,-) dg{x,y) dh{x,y) 

+ sup ^^^-^ sup sup sup 

Od-hiki,-) ^a-gi-M) dy dx 
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+[sup — Y sup ^ 

ox da-g(-M) 

> sup — Ysup- — ^^ + 2 sup — sup — 

oy Od-hiku-) oy 

dh{x,y) ^ 9d-h(o,-) 
+ [sup — ] sup- ^ 

dX ^a-g{-M) 

since 9a-g{.fi) > Oa-g{;k2)^^d-Ho,-) > Od-h{ku-)- Therefore, by the Theorem 3.1, 
(13) has a unique coexistence state for aU (a, rf) G dV . Furthermore, by the 
estimate of the solution in the proof of Theorem 3.1, if {u,v) is a positive 
solution for (a, rf) G F, then 

.^^^ dh{x,y) 
B dx B dy 
^ r dg{x,y) dh{x,y) 9d-h(o,.) dg{x,y) Oa-g(.fi) 

> [sup + sup sup [sup sup 



+ sup 



dy dx Oa-g{.,k2) dy 9d-h{ku- 

dh{x,y) 
dx ^ 



^ r dg(x,y) dh{x,y)v dg(x,y)u dh{x,y) 

> sup — h sup — sup — + sup — . 

oy ox u oy V ox 

Thus, we obtain 

dg{x,y) . dh{x,y) dg{x,y) dh{x,y) ^ 

4 mf — mf — - — -uv > sup — u + sup — — -v\ . 

B Ox B Oy Oy Ox 

It implies that if {u,v) is a positive solution of (13) for {a,d) G F, then the 
Prechet derivative of (13) at {u, v) is invertible from Theorem 4.3. Therefore, 
the theorem follows from Theorem 5.1. 
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